clc;
clear;

V = [0, 0; 1, 0; 0, 1];

[gauss_points_loc, gauss_weights_loc] = gauss_lengdre_loc(V, "2d");

phi = basis_function_loc(V, gauss_points_loc(:,1), gauss_points_loc(:,1), "P2", "dx");

[phi{1}, phi{2}, phi{3}, phi{4}, phi{5}, phi{6}]